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Abstract 

We present an improved partial wave analysis of the dominant LSP annihilation 
channel to a fermion-antifermion pair which avoids the non-relativistic expansion 
being therefore applicable near thresholds and poles. The method we develop allows 
of contributions of any partial wave in the total angular momentum J in contrast 
to partial wave analyses in terms of the orbital angular momentum L of the initial 
state, which is usually truncated to p-waves, and yields very accurate results. The 
method is formulated in such a way as to allow easy handling of CP- violating 
phases residing in supersymmetric parameters. We apply this refined partial wave 
technique in order to calculate the neutralino relic density in the constrained MSSM 
(CMSSM) in the presence of CP-violating terms occurring in the Higgs - mixing 
parameter /i and trilinear A coupling for large tan (5. The inclusion of CP-violating 
phases in /x and A does not upset significantly the picture and the annihilation 
of the LSP's to a bb, through Higgs exchange, is still the dominant mechanism in 
obtaining cosmologically acceptable neutralino relic densities in regions far from the 
stau-coannihilation and the 'focus point'. Significant changes can occur if we allow 
for phases in the gaugino masses and in particular the gluino mass. 



1 Introduction 



In the framework of Supersymmetric Theories, in which R - parity is conserved, the 
lightest of the neutralinos, x, which is also the lightest supersymmetric particle (LSP) is 
a promising candidate for cold Dark Matter (DM) [1]. 

The knowledge of its relic density today combined with recent experimental observa- 
tions, especially those coming from WMAP [2,3], put severe constraints on the parameters 
of the theory which if combined with other data enhance the possibilities of discovering 
Supersymmetry in future accelerator experiments. 

In the framework of the Constrained Minimal Supersymmetric Standard Model small 
values of the neutralino relic density, m the vicinity of 10% range as recent as- 

trophysical data on Cold Dark Matter suggest [4-6], [3], are obtained in regions of the 
parameter space which belong to either the neutralino - stau coannihilation region [7, 8] , 
or to the focus point [9, 10], or to a region in the neighborhood of the line along which 
the pseudoscalar Higgs has a mass approximately twice that of the LSP mass [5, 11-16], 
which is characterized by large tan/?. In this region the LSP annihilation fusion to a 
pseudoscalar Higgs, which subsequently decays into a bb pair, dominates the annihilation 
process. Interestingly enough in the same large tan/3 region the elastic LSP - nucleon 
cross section receives one of its highest possible values close to the sensitivity limits put 
by future experiments which will directly search for DM [15, 17]. In this region as well 
as in the neutralino - stau coannihilation region the LSP is mainly bino. In the focus 
point region the neutralino LSP bears a relatively sizeable higgsino component and its 
pair annihilation cross section through the Z and/or the light Higgs boson exchange be- 
comes significant, resulting to cosmologically acceptable relic density and high values for 
the LSP-nucleon cross section. 

In order to calculate the LSP relic density, one has to solve the Boltzmann equation 
having as input the thermal average (crv re i) of the LSP's annihilation cross section [1, 
18-30] cr(s) times their relative velocity v re \ , if we are in a region of the parameter 
space where coannihilation effects are negligible. Its computation demands integration 
over the center of mass scattering angle cos Oqm of the transition matrix element squared 
summed over the final spins followed by an additional integration over s through which 
the aforementioned thermal average is defined. 

If one follows the standard trace technique for calculating the amplitudes many in- 
terference terms are encountered since the squares of the annihilation amplitudes are 
summed in a coherent way. However this can be successfully accomplished as has been 
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shown elsewhere [28-30]. 

In this paper we shall pursue a partial wave analysis in the total angular momentum 
J , which can be considered as complementary to the trace technique and it improves 
analogous treatments existing in literature. In a partial wave expansion the integration 
over the scattering angle cos 9qm is carried out trivially, since the amplitudes of different 
J are summed incoherently in the total cross section avoiding numerous interference 
terms. This allows a better theoretical and numerical control of the annihilation process. 

In addition the analysis , which we employ in this paper, has the following advantages: 

i) It avoids the non-relativistic expansion resulting to more accurate numerical treat- 
ment of the thermally averaged cross section especially near poles some of which, as 
those of Higgs bosons, are significant in obtaining small LSP relic densities in accord 
with recent cosmological data. The partial wave schemes encountered in literature 
usually expand the total cross section in powers of the relative velocity v re i keeping 
terms up to 0(v^ el ) 1 . 

ii) It allows contributions of partial waves of any J, where J is the total angular 
momentum, offering a better approximation scheme of the partial wave series. The 
existing partial wave treatments are usually truncated to include terms up to p- 
waves in the orbital angular momentum L of the initial state which are included 
in J < 2 . 

iii) The scheme is formulated in such a way as to allow easy handling of CP- violations 
residing in the superpotential Higgs mixing parameter /x and the trilinear A soft 
couplings and/or in the presence of additional supersymmetric CP-violating phases 
occurring in nonminimal models. 

The scheme we present can be directly compared to the exact results one obtains 

using the trace technique and from this the importance of terms larger than J = 2, 

omitted in other approaches, can be sought. As a preview, our analysis reveals that 

the thermal average of the cross sections times the relative velocity of the annihilating 

LSP's obtained using J < 2 in the partial wave series approximates the exact result to 

an astonishingly good accuracy in the entire region of the constrained MSSM parameter 

space for temperatures that are relevant for the calculation of the relic density. We shall 

comment on this later on. 

1 In the first of reference [11] such an expansion holds and the problems associated with the poles and 
thresholds are avoided by not expanding the kinematical factors /?/ and the s-channel propagators. 
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In this paper paying special emphasis to the mechanism for the LSP annihilation 
through the Higgs boson exchange, we apply this refined technique in order to calculate 
the neutralino relic density in the CMSSM in the presence of CP-violating terms residing 
in the \x and trilinear A parameters of the CMSSM [31] taking also into account the effect 
of the mixing occurring in the Higgs sector. The effect of the coannihilation processes, 
in regions of the parameter space where this applies, in the context of this partial wave 
approach, will be the subject of a forthcoming publication. 

This paper is organized as follows: In section 2 we give a brief outline of the helicity 
technique. In subsection 2.1 we discuss the s- channel contributions of Higgs and Z - 
boson exchanges while in 2.2 we present an analysis of the t and u - channel sfermion 
exchanges. In section 3 we briefly discuss the handling of the thermally averaged cross 
section and in section 4 we apply this technique when CP-violating phases are present. 
We end up with the conclusions presented in section 5. 



2 The Helicity Formalism 



In this section we will briefly outline the helicity formalism relevant to our calculation 
of the amplitude for the annihilation of two LSP's to a fermion - antifermion pair. We 
should evaluate the annihilation amplitudes for the above processes distinguishing initial 
and final helicity states . The Dirac spinor of a free fermion with mass m , 3-momentum 
p and energy E is: 



u{p) = N 



E + m 



(i) 



where £ is a two-component spinor. The overall normalization factor N = 



,E 



m 



2m 



V2 : 



is 



consistent with uu—1. In order to find the spin - \ helicity states , we choose the two 
component spinor £ to be eigenstates of the helicity operator A = \ a ■ p: 

1 



A 6 = A6 



(2) 



p is the unit momentum vector, p 
states are found to be, 
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corresponding to spinors of helicities +~ and — ~ respectively. Therefore the one - particle 
helicity states for a particle, u(p, ±), or antiparticle, v(p, ±), with helicities ±| are given 
by 



u{p, ±) 



N 



IpI 



E + m 



± 



v(p, ±) = N ( E + m 



6 



(4) 



Note that u(p, ±) = Cu T (p, ±). 

For the construction of a two - particle helicity state, one has to define appropriately 
the helicity spinors for the second fermion, i.e for the fermion carrying momentum — p 
in the center of mass frame. If the momentum of the first particle p is pointing in the 
direction specified by the angles (9, <p), then the momentum of the second particle is in the 
direction specified by (it — 9, (p + n). The corresponding helicity spinors for the second 
fermion are, 



/ 



V 



-tip 



sin 



cos- 



and 
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cos- 
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(5) 



Sm 2 / 



To comply with the Jacob- Wick phase convention these have been multiplied by the 
appropriate phase factors — 2\e~ l2Xlp 2 . 

Any helicity amplitude can then be expressed as a sum over the total angular momen- 
tum J, 

•MasA^M,^) =E( 2J+1 )4>) MaUa^K ^ ( 6 ) 



is the reduced matrix 



where \i = Ai — A2 and // = A3 — A4. In this equation M. 
element, and the rotational functions d^,(9) satisfy the orthogonality relations 



(J) 

A3A4;AiA2 



+ 1 



d(cos9)d^;(9)d%}(9) 



2J 



(7) 



so the interference terms are avoided when integrating the amplitudes squared. 

We are primarily interested in studying the annihilation of two neutralinos into a 
fermion-antifermion pair using the helicity formalism. The relevant channels concern 
Z, Higgs (s-channel) and sfermion / (t and u channel) exchange, see figure d For every 
helicity amplitude we need calculate the following matrix elements involving the spinors 
of the initial, X\i (Pi) > Xx 2 (P2), an d final states, (p 3 ) ,/a 4 (P4), which are denoted by 
2 For more details see for instance [32] and references therein. 
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<0| r |xa 1 (Pi)xa 2 (P2)> , (f\ a (P3)f\ 4 (P4)\ r |0) , (A 3 (P3)|r|x Al2 (pi )2 )), and 
(fxiiPi) I r Xai, 2(^1,2)) • These are the helicity dependent fundamental building blocks 
encountered in the vertices of the s, t and u— channel Feynman diagrams. In these V 
stands for Y = {/, 7 M , 75, 7^75} depending on the case under consideration. For lack of 
space we do not present their analytic expressions. 

Each of these blocks is expressed in terms of the helicity spinors described previously. 
These are functions of the angle 6, the scattering angle in the center of mass frame, 
assuming that the two annihilated neutralinos move along the z— axis, the azimuthal angle 
and the center of mass energy squared s. In the center of mass frame the magnitudes 
of the incoming and outgoing particle 3-momenta are, 

PiM = (--m ji ) , p 3 , 4 (s) = (- - m f ) . 
and the Mandelstam variables t and u as functions of s and 9 are given by, 

t = m| + mj - I + 2px(s)p 3 (s) cos# 

OrS /\/\ r\ 

u = + rrij- — — — 2pi(s)p 3 (s) cosv . 

In this formalism the unpolarized total cross section for the annihilation of the two neu- 
tralinos into a fermion - antifermion pair, in the center of mass frame, is given by 

"W = 6^T7E E pj+dKUj' (8) 

J A3A4;AiA2 

and hence no integration over the scattering angle is required once the helicity amplitudes 
M-x^x A -x \ 2 are known. Loss of accuracy occurs by truncating the partial wave series 
to a maximum value J m ax of the total angular momentum. However as stated in the 
introduction the series converges fast in the regime of energies which are relevant for the 
calculation of the relic density. In this paper we use a value of J ma x — 3 which is good 
enough to obtain accurate results and at the same time a fast code which returns the 
cross section and its thermal average. 

Using the machinery outlined before one can express any s, t and u— channel scattering 
amplitude as functions of s and cos# and bring them in the form given by Eq. (jHJ. This 
will provide us with the helicity amplitudes appearing in Eq. (JHJ). The analytic expressions 
for the s, t and u- channel helicity amplitudes are listed below. 

2.1 s-channel contributions 

The s - channel are relatively easy to handle. Since the exchanged particle carries spin of 
either s — 1, for the Z - boson, or s = for the Higgses, the partial wave series terminates 
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Figure 1: Graphs that contribute to the process XX —> f f- Hi in the second graph denote 
the Higgs mass eigenstates 

at J < 1. 



2.1.1 Z-exchange 

The reduced matrix element for the Z - exchange, as it follows from the Feynman rules 
using the helicity blocks described in the previous section, is 



M (J),z 

A3 -Vi;AiA2 



s - M 2 z + iM z T z 



KxP^ (V f FU + AfF xk) A xx 



(9) 



where the vertices N x x , F XaXi and F x ^ Xi can be expressed in terms of the following matrix 
elements, 



<A 2 
A3A4 



(0|7 M 7s I Xx 1 (Pi)X\ 2 (P2)} 
(fx s (ps)fxM) It* |o> 
(fxMfxM |7 M 75 I 0) 



A xx is the coupling of the Z— boson to the neutralinos and V*, A*, are the vector and 
axial couplings of the Z— boson to a fermion - antifermion pair labeled by /. 

k k 

The numerator of the Z propagator, in the unitary gauge, is P^ u = —g^ u + p ' £ where 

M z 

9fiv — diag(l, —1, —1, —1). Since the four-momentum of the exchanged Z- boson in the 
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center of mass frame is = (y/s, 0) the tensor P^ u receives the following form 



p 



~ 1 + W Z ' fi = v = 

, n — 0, v — % 



~9ij , fi = t,v = j . 
where i, j = 1, 2, 3. 

After evaluating (jUJ) , we find that a non-zero contribution emerges only when the 
helicities of initial and final states are combined to give J = or 1 as expected. In 
fact 



^x 3 x,MX 2 ~ s -Ml + iM z T z K ^-MX 2 - 



with 

= (-l) A3 " Al ^oVo(4m,m / )(-l + ^)^ 

z 



Vo (2m/) y/2 V f + 2 (<Vi + fy-i) (^ ^ + A/ s ~ 4m f ^) 



111 



2.1.2 Higgs exchange 

Using the technique outlined previously one can obtain the corresponding amplitudes for 
the s - channel annihilation process through the Higgs mass eigenstates H iy i = 1,2,3 

M XaXi . MX2 - S -M 2 Hi + iM H T Hi A ^A 4 ;A 1 A 2 W 

In this equation , the superscript Hi denotes the exchanged Higgs and the non - vanishing 

(J) H 

■^x A4 A1A2 are those corresponding to J = only, since Higgses carry zero spin. These are 
given by 

R (J=0), Hi = 
A3A4;AiA2 

<WW (#^-4m| + ^v^(-l) Al ^) (^y s -4m 2 / + ^yi(-l) A3 ^) (13) 

In these , g^^g^ denote the scalar(S') and the pseudoscalar (A) couplings of the neu- 
tralino LSP x to the Higgs H i: while g$ A are the corresponding couplings of the same 
Higss to a fermion pair //, 

~ Xigf + gf 7 5 )X H { + f(g f / + g f / 75)/ H* 
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Note that in our final results we should include a minus sign in the s - channel amplitudes 
to account for the relative (-1) difference from the corresponding t-channel amplitudes. 

Note that we have allowed for the more general coupling of Higgses to fermions and 
the lightest of the neutralinos so that our formulae are applicable when supersymmetric 
CP- violating phases are present in the A and fi parameters which in turn cause mixing 
of the CP-even and CP-odd Higgs eigenstates. 

2.2 t &; u-channel : sfermions exchange 

The t and u amplitudes are rather difficult to evaluate since the exchanged sfermion 
propagators depend on the scattering angle 6. The propagators can be expanded using 
the formula 

— !— = V U n {a)P n {x) (14) 
1 — ax z — ' 

where P n (x) is the first Legendre function of order n and its argument is x = cos#. In 
this equation a < 1 and 

tt r\ 2n + l 7 P w (cosfl) 

n « 0) = / ~, -d(cos0). 15 

2 J 1 — a cos v 

-i 

The function H n (a) appearing in the above formula is actually (2n + 1) Q n (l/a)/a where 
Q n is the second Legendre function. Multiplying by the corresponding vertex contribu- 
tions, the t and u channels can be cast in the form (jOJ). Instead of presenting separately 
the contributions of the t and u channels we find it more convenient to present their sum. 
Note that a relative (-1) sign between the t and u channels must be included which has 
been duly taken into account in the following expressions. We find that their sum gives 
rise to the following partial wave amplitudes 

u(J)J _ \- (2J+ 1)' 1 / , , x f g l s * \ nfi x 

- /W A 3 A 4 ;A 1 A 2 - 2^ F^S) ^ J^A 3 A 4 ;A 1 A 2 ^ _ g i* J l iD J 

/j. The constants g l si9\ are the scalar and pseudoscalar neutralino-fermion-sfermion 
couplings defined as 

M + 9 l Al 5 )f!f + h.c. 



while the quantity Fi(s) is given by 



2,2 2 

m~ + m f — m? 

X J Ji 
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The matrices -^a 3 a 4 a x a 2 depend on the energy E = \fs/2 through 



A± = \ (E + m f )(E + mz) ±J(E-m f )(E- m 2 ) 



B± = \ (E + m f )(E - mz) ± J (E - m f ) (E + 



m 



X) 



2pm y (s - Am 2 f )(s - 4m 2 -) 

and also on the functions HjiaA whose arguments are a* = - . . = — t——. 

S Fi(s) 2F i {s) 

For large energies tij approaches unity and ITj approaches a logarithmic singularity. In 
table £Q the matrices F^^. XiX2 are presented for all helicity combinations of the incoming 
LSP's and the outgoing fermions. 



3 Handling the thermal average and the relic density 

The primarily task in calculating the relic density is to find the thermal integral 3 

1 If 00 

(™ re l) (T) = — y- — - — — - / dspWtfKxiy/i/T). (17) 

The momentum p is related to the CM energy and through s = 4 (p 2 + m|). Scaling 
the temperature in terms of we can use as thermal variable x = ^7. Introducing a 
new variable y defined by 

9 =i(^-2), (18) 

Eq. ()17|) can be cast in the form 

1 1 f°° 

(av rel ) (x) = / dy (xy + 2) ^(^ + 4) W(y) K^y + 2/x) . (19) 

In this equation is defined to be W(s) with s replaced by its value from Eq. 

that is v^i = m^(xy + 2). Before continuing we should perhaps comment on how the 
singular behavior of the Bessel functions encountered in (|17jl. or same (|19jl. are evaded 
in our numerical code. This problem is not related to the partial wave expansion we 
employ to calculate the cross sections but it is due solely to the behavior of the Bessel 
functions appearing outside and within the integrands in the above relations. In fact the 
presence of the /^(l/x) in the denominator of Eq. (fTTj) is potentially dangerous, since 
for x < x c , with x c a small number of order 0.003 or so, the Bessel function underflows 
3 The W(s) that used here is related to the cross section by a(s) = W(s)/py/s (see Ref. [19]). 
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K 2 ~ 0, resulting to numerical overflow. However, the asymptotic expansion of the Bessel 
functions for large arguments is known to be 



K n (z)^J^e- z P n (z), (20) 



where c = l/x c and 



( te»-i)(4 B '-a') 

* 1182 2!(8z) 2 + "' 1 ' 

Using Eq. (}2T)j) one can obtain a form of Eq. (fT§|) which is free of such overflows in the 
offending region x < x c 



1/1 I f 00 

(av rel )(x) — > — W— / cfy v/y(xy + 2)(xy + 4) W(y) e"* P^y + 2/x) . 

(22) 

For x > x c the K% does not underflow but the Bessel function K\ within the integral does 
in the integration region in which its argument exceeds l/x c . In this integration region 
and in order to avoid such a numerical underflow we approximate K\ by Eq. (J2*U|) within 
the integral. These results can be applied for an accurate numerical evaluation of the 
integral of Eq. (|19|) valid for any x. 

In obtaining the LSP relic density one has to solve the Boltzmann transport equation 
numerically or use an approximate solution which however is quite accurate. Regarding 
the first approach the details of how one solves the Boltzmann equation are described 
by Edsjo and Gondolo and by Gelmini and Gondolo in reference [27]. Also an efficient 
method, reminiscent of the WKB approximation, for solving Boltzmann equation is pre- 
sented in [23]. The details concerning the second approach, which is mostly used in the 
literature, are given in [1,33], or by Griest et al. in the second of references in [21]. Ac- 
cording to it the freeze - out temperature Tf and the corresponding point Xf is determined 
by 

\ -i -i 

0.038 g eff Mpi anck (av re i) 



In 



1/2 




(23) 



i li 

In this equation g e ff are the effective energy degrees of freedom and gj is the total 
number of effective relativistic degrees of freedom related to the entropy degrees of freedom 
h and g eff by 

1/9 ( , dh \ 1 i , 

°- /2 = { h+ irjr)^jj- < 24 > 
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Figure 2: The ratio ((av re i) 3 — (av re i) 2 ) / (av re i) 2 , in %, as a function of x. The 
subscripts refer to the value of J max in each case. The figure corresponds to values 
A = 0, m = 1.2 TeV, tan/3 = 55 and values of mi/ 2 equal to 900 GeV (solid line) and 
800 GeV (dashed line). 

Eq. (J23)) can be iteratively solved to yield the freeze - out value Xf and from this by 
integrating (av re i) (x) one can have the integral 



The value of the freeze - out point Xf is roughly ^1/20 and thus the values of (av re i) 
needed for for the calculation of (f2*U|) refer to points x < Xf as is obvious from (J2HJ). 

In order to estimate the error induced by truncating the partial wave series to values 
J < Jmax in figure El we display the ratio ({<Tv re i) 3 — (<?v re i) 2 )/ (av re i) 2 , where the 
subscripts refer to the value of J ma x in each case. The plot of the figure corresponds to 
values A = 0, m = 1.2 TeV, tan/? = 55 and values of mi/ 2 equal to 900 GeV (solid 
line) and 800 GeV (dashed line). Both cases refer to fi > 0. The figure clearly shows 
that this ratio is small ~ 10 -4 when x lies between and 0.5, which includes the range 
of integration in Eq. |23 indicating that the partial wave series saturates the exact result 
if one truncates at a value J max = 2. This is a general feature valid in the entire region of 




(25) 



The LSP relic density is then given by the expression 




(26) 
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the parameter space. In our numerical analyses concerning this work the value of J max-, 
which is actually input in our code, is set to J max — 3. Higher values will yield more 
precise results at the expense of slowing down our numerical code. 

It should be pointed out that the partial wave expansion we employ in terms of the 
total angular momentum J differs from the expansion in terms of the orbital angular 
momentum L. Since L + S is even for the initial state of the two Majorana LSPs, 
expanding up to J = 2 includes not only the L = 1 terms (p-waves) but in addition all 
amplitudes characterized by L = 2. This is much improvement over the existing orbital 
angular momentum schemes which are truncated to L — 1 equivalent to keeping terms 
up to 0(v^ el ) in the relative velocity of the two LSPs. For this reason the series in terms 
of J yields very accurate results if already truncated at J max = 2, as shown in figure El 

4 Neutralino abundance with CP- violation at large 

tan (3 

We apply the above described helicity amplitude technique to calculate the neutralino relic 
density at the large tan/3 regime, where it is known that in the CP-conserving case the 
contribution for the neutralino pair annihilation cross-section through the pseudo-scalar 
Higgs boson A exchange plays an important role in obtaining cosmologically acceptable 
relic densities far from the coannihilation and focus point regions. Since in this work we 
are mainly interested in the aforementioned mechanism for the annihilation of the LSP's, 
other production channels, e.g. gauge, Higgs bosons etc, do not contribute significantly. 
Hence we consider the neutralino LSP pair annihilation to fermion pairs at the relativistic 
level, including all finite width contributions. The code we employ, which is based on the 
partial wave analysis we outlined in the previous chapters, is pretty fast and can cope with 
the Higgs resonances which play an important role in obtaining an LSP relic abundance 
in accord with the recent cosmological data. As already stated in the previous section the 
value of the maximum allowed angular momentum J max is input but a value of J max = 3 
is sufficient in obtaining a very good accuracy. 

In this work we have focused our interest in the so-called "funnel" region of the pa- 
rameter space, where one obtains cosmological accepted relic density through annihilation 
to Higgs boson resonances. In principle one expects that relatively sizable CP- violating 
effect can also occur in the focus point region, where the LSP carries a non-zero Higgsino 
component. This results to an enhancement of the LSP pair annihilation cross section to 
massless fermions, through the Z-boson exchange. Although we will not present numeri- 
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m 1/2 (GeV) m 1/2 (GeV) 

Figure 3: The neutralino relic density, in the (mo,mi/2) plane, for tan/5 = 50 and 55. 

cal results in the focus point region, we have checked that there the CP- violating effects 
are almost comparable with those around the Higgs pole. Furthermore, the latest CDF 
and DO data [34] suggest a somewhat higher value for the top quark mass m t = 178 GeV, 
pushing the focus point region up to considerably high values of m . 

The method we develop is capable of accommodating CP-violating sources residing 
in the \i and the trilinear A parameters as well as in other parameters, as for instance 
in the gaugino mass parameters, which are important when we depart from the minimal 
schemes. In this work we consider the constrained scenario and hence only the fi and A 
phases are allowed. In figure EH and for nonvanishing phases equal to </> M , 4>a = tt/3 at the 
weak scale we present in the mo,m 1( / 2 plane the allowed domains (in green), which are 
consistent with the recent cosmological data from WMAP [2] % hi = OA126tomt\ • The 
cases shown are for values of tan /3 equal to 50 and 55 respectively. The shaded area (in 
magenta) at the bottom, in both panels, is excluded since in this the stau is the LSP. In 
figure 0] we display the relic density for fixed mo, A Q as a function of for tan/3 = 55 
and for nonvanishing phases <f)^,(f)A — tt/3 of fi and A. For comparison the case with 
<Pii) 'Pa — is also shown. In the left and right panel of this figure the cases A = and 
Aq = —500 GeV are presented. We should remark that in our analysis we have relaxed 
the EDM constraints [35, 36] in order to locate possible regions of the parameter space 
which yield substantial deviations from the results of the CP-conserving case. 

From figure |U one sees that no significant change is observed and the results for the 
CP-conserving and CP-violating cases are almost identical. This holds for other values 
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£2 Z h 2 tanp=55, A = 0, m =1200 GeV, n>0 Cl x h 2 , tanp=55, A = -500 GeV, m =1200 GeV, |i>0 




500 750 1000 1250 1500 500 750 1000 1250 1500 

m 1/2 (GeV) m v2 (GeV) 

Figure 4: The relic density for fixed mo, as a function of mi/2 for values of the 

phases and the remaining parameters shown in the figure. Only the /x > cases are 
shown. For comparison the corresponding values when the phases are zero are displayed 
(solid line). The horizontal dotted lines mark the cosmologically allowed region from the 
WMAP data. 

of m , A and phases as well showing that the effect of the phases plays no important 
role to the relic density. Phrased in another way, despite the fact that the CP-even and 
CP-odd Higgs eigenstates get mixed when CP is violated, the Higgs exchange is still the 
dominant mechanism in obtaining small values of the relic density, as long as we are away 
from the coannihilation and focus point regions. The mixing of Higgses [37, 38] induced 
by the phases of fi and A has little effect and this is in agreement with the findings of 
Ref. [39]. Thus the CP-conserving case results are not destabilized by switching on the 
phases of M , (pA parameters and the conclusions reached in previous analyses [5,12,15,16] 
hold true. 

Significant changes in the relic density can occur if one allows for the existence of 
phases in the gaugino mass parameters mainly because the latter affect the value of the 
bottom quark mass as shown in [39]. In particular, the supersymmetric QCD corrections 
to the bottom quark mass depend on the phase £3 of the gluino mass M 3 at the GUT scale. 
Even if we ignore the phases M , <pA as well as the phases of the gaugino masses M 1; M 2 
and the angle specifying the misalignment of the vev's of the Higgs fields, the effect of £3 
is quite important and can drastically change the picture. In the discussion that follows 
we have only allowed for a nonvanishing £3 setting all other phases to zero. With only 
£3 present the electroweak sector is not directly affected by it, up to the one loop order. 
At this order the phase £3 affects the electroweak sector implicitly through changes in 
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the bottom Yukawa coupling as described earlier. In particular the Higgs sector is the 
same as in the CP- conserving case and the 1-loop expressions for Higgs masses still hold. 
It is known that the mass of the pseudoscalar Higgs A is quite sensitive to the size of 
low energy threshold supersymmetric corrections to the bottom Yukawa coupling [16,39]. 
Consequently £3 influences the position of the A-Higgs pole, M4 ~ 2m^, where the rapid 
neutralino pair annihilation yields cosmologically acceptable relic densities. In order to 
see the effect of the phase of M 3 on the relic density, in fig. El we plot the values of the 
pseudoscalar mass Ma , the value of the mass of the annihilated neutralino pair 2m^ and 
the relic density as functions of the phase £3 for two characteristic sets of input parameters 
shown on the left and right panel. As already stated we have only switched on the phase 
£3. On the left panel it is observed that although the value of tan/? is not quite high 
the approach to the pole, Ma — 2m^, takes place for values £3 ~ 37r/4 resulting to rapid 
X annihilation through A-Higgs boson exchange and small relic densities in accord the 
WMAP data ^0 — 0.1). On the right panel and for different input values we observe 
that the relic density becomes significantly smaller for values £ ~ n, but not compatible 
with the WMAP bound. In both case the plots look symmetric upon £ 3 — > — £3. This is 
due to the fact that the SQCD corrections are insensitive to the sign of £3 if the other 
phases are put to zero. 

In the left panel of figure El we plot the neutralino relic density for various values 
of £ 3 and the same input parameters m = 1200 GeV, A = 0, tan (3 = 55, // > 0. One 
can see that for £ 3 = 0, CP-conserving case, the minimum of the fi,% appears at 
mi/2 = 1000 GeV, while for £3 = tt/G this occurs at a lower value m\j% = 750 GeV. This 
means that as £3 increases, the funnel of the right panel of figure El is widening up and 
moving to the left allowing for smaller values of m-1/2 and hence lighter supersymmetric 
mass spectrum in general. This is clearly seen in the right panel of figure El which should 
be compared to the right panel of figure El 4 - 

5 Conclusions and Outlook 

We have presented an improved partial wave expansion technique for calculating the neu- 
tralino pair annihilation to a fermion pair, based on the helicity amplitude method. In 
the partial wave treatment the amplitudes carrying different total angular momenta are 

4 In figure l^the phases <j)^, cf>A have been taken equal to 7r/3 and not vanishing as in figure HO However 
when £3 is zero the effect of (j)^, 4>a is negligible and the figure [31 is the same with the corresponding one 
in which t/^, (j>A are put to zero. 
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Figure 5: The mass Ma (solid line), (dashed line) and the neutralino relic density 
(dotted line) as functions of £3 when all other phases are set to zero. In the left (right) 
panel m = = 800 GeV, tan/3 = 35 (m = = 500 GeV, tan/5 = 30), and 

A = 0,/i > 0. 

summed incoherently in the total cross section, avoiding numerous interference terms, 
resulting to a better theoretical and numerical control of the neutralino annihilation pro- 
cess. Also the non relativistic expansion is avoided obtaining a more accurate numerical 
treatment of the thermally averaged cross section especially near poles some of which, as 
those of Higgses, are significant in obtaining small LSP relic densities as recent cosmo- 
logical data suggest. The scheme we employ allows for contributions of partial waves of 
any J, offering a better approximation of the partial wave series. The methods existing 
in the literature are usually truncated to include terms up to p - waves in the orbital 
angular momentum L of the initial state of the annihilating LSP's corresponding to val- 
ues of J < 2. Our scheme is formulated in such a way as to allow easy handling of 
CP-violations residing in the superpotential Higgs mixing parameter [i and the trilinear 
A soft couplings as well as in the presence of additional CP- violating phases occurring in 
nonminimal extensions of the MSSM. 

Using this method we include non - vanishing CP-violating phases M , </>a to the /z and 
A parameters of the CMSSM. We study the effect of these phases at the large tan (3 regime, 
where it is known that the exchanges of Higgses is one of the dominant mechanisms to 
reduce the values of the LSP relic density to cosmologically acceptable levels if CP is 
conserved. We found that even if one relaxes the EDM constraints the effect of these 
phases is small, in the entire region of the parameter space away from the coannihilation 
and focus point regions, and the values of the neutralino relic density obtained are almost 
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Figure 6: (Left panel) The neutralino relic density as a function of m^, for m = 
1200 GeV, A = 0, tan/3 = 55,/i > 0, for £ 3 = 0, n/8 and 7r/6 (solid, dashed and dashed- 
dotted lines). The WMAP favored region is marked by the horizontal dotted lines. 
(Right panel) The neutralino relic density, in the (mo,m^ 2 ) plane, for tan/? = 55 and £ 3 
equal to x/6. 

the same as in the CP-conserving case. Therefore the Higgs exchange mechanism is still 
one of the dominant mechanism in obtaining small values for the relic density and the 
CP-conserving case results are not upset by switching on the phases (p^, 4>a- Significant 
changes in the relic density can occur if one allows for the existence of phases in the gaugino 
mass parameters as shown elsewhere [39,40] which can only happen if one departs from 
the minimal scenario. 

The method of partial waves presented here will be generalized to cover the subdom- 
inant processes of annihilation of the neutralino LSP's to channels other than that of a 
fermion pair and to include also the coannihilation processes which are important in par- 
ticular regions of the parameter space. The results of such an analysis will be presented 
in a forthcoming publication. 
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Table 1 



A3A4A1A2 



Mi 



A3A4;AiA2 



+ + + + 

+ + 

- - + + 



fx' = 



M 



(A3A4) 

Ai 



J even 



n J + (-i) A 3-Ai(n J _ 1 



2J- 1 



n J+ i 



J + i , 

2J + 3' 



J even 



- + ++ 



// = 



(A3A4) 



\M</ + i)( 



n 



n 



j+i 



2J-1 2J + 3' 



J > 1 



/^0,// = 



i g 2 (A 3 A 4 ) + (.^Jj^ 



(A3A4) 



n 



j-i 



n 



j+i 



2J-1 2J + 3' 



(— 1) J eg) (previous) 



+- 



i ^ 2 (A 3 A 4 ) + (_ 1 )J^A 3 A 4 ) ^ 



J > 1 

(A3A4) , 1 



-1) J M 3 
J + l 



(A3A4) 



n,+ 



LJ-1 



2J- 1 



n 



.7 



j+i 



2J + 3' 



-1) J <S> (previous) 



Table 1: The matrices F appearing in Eq. pfij) . The matrices M are given in table El 
The ± in the helicity column stands for ±=. 
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Table 2 



H = 0,// = 



M 



++ 



Ml 



Ml 



M 

-A_B_ 

-A + B + 

A + B + 
-A 2 _ 
A_B_ 



-A_B_ 

-Bl 
-A + B + 

A + B, " 

A_B_ 
-B 2 



M: 



Ml 



+- 



Afr + 



-B_B + 
A + B_ 
-B_B+ 
A + B_ 
B_B + 
-A_B 

B_B: 



Ml 



Ml 



Ml 



Ml 



A^A + 
A + B_ 

A_B+ 

A_B + 

A,iL 



-A_B + 
-B_B + 
-A+ZL 

-A+lL 
B_B + 

-A_£? 4 
B_B+ 



M+- 



M 3 + " 



-A 2 




-A_iL 






-A+B+ 








A + B + 


-A+5+ 


B\ 


-A 2 _ 


-A_B_ 


A_zL 


B 2 



Table 2: The matrices M appearing in tabled 
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